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ABSTRACT 

The number of strong (equivalent width > 1A) Mgll absorbers observed towards Gamma-ray bursts 
(GRBs) has been found to be statistically larger than the number of strong absorbers towards quasi- 
stellar objects (QSOs). We formalise this "Mgll problem" and present a detailed explanation of the 
statistical tools required to assess the significance of the discrepancy. We find the problem exists at 
the 4cr level for GRBs with high-resolution spectra. It has been suggested that the discrepancy can be 
resolved by the combination of a dust obscuration bias towards QSOs, and a strong gravitational lens- 
ing bias towards GRBs. We investiga te one of the tw o most probable lensed GRBs that we presented 



in our previous work (GRB020405; | Rapoport et al. ) and find it not to be strongly gravitationally 
lensed, constraining the percentage of lensed GKBs to be < 35% (2a). Dust obscuration of QSOs has 
been estimated to be a significant effect with dusty Mgll systems removing ^20% of absorbed objects 
from flux-limited QSO samples. We find that if ~30% of the strong Mgll systems towards QSOs are 
missing from the observed samples, then GRBs and QSOs would have comparable numbers of ab- 
sorbers per unit redshift. Thus, gravitational lensing bias is likely to make only a modest contribution 
to solving the Mgll problem. However, if the dust obscuration bias has been slightly underestimated, 
the Mgll problem would no longer persist. 

Subject headings: Gamma-ray bursts: general, Gamma-ray burst: individual: GRB020405, Gravita- 
tional lensing: strong 



1. INTRODUCTION 

The number of strong Mgll absorbers per unit red- 
shift along Gamma-ray burst (GRB) lines-of-sight (LOS) 
has been found to be statistically differe nt than the rate 
towards quasi-stellar objects (QSOs). |Prochter et al.| 
(2006bJ hereafter P06) fo und 14 strong absorbers along 
the LOS to 14 GRBs, and |Vergani et al | fl2009| hereafter 
V09) increased the sample and found 22 strong absorbers 
along the LOS to 26 GRBs. In contrast, studies of QSOs 
have found Mgll absor bers in only ~ 25% of sight lines 
(Prochter et al. 2006a). Both being high-redshift bea- 
cons, GRBs and QSOs are expected to have similar LOS 
through the cosmos, and explaining the preponderance 
of Mgll absorbers towards GRBs has pro ven a challenge. 

In t he effort to unveil this problem, Porci anT et al.| 
2007| ) considered dust obscuration, beam size differ- 
ences, intrinsic properties of GRBs and gravitational 
lensing as possible causes. They found beam size dif- 
ferences irrelevant, with simulations predicting the ab- 
sorbing systems are sign ificantly larger than either beam. 
Cucchiara et al.| ( |2009[ ) found no dissimilarities between 
the populations ol absorbers, suggesting it is unlikely 
that there are excess absorbers physically associated with 
the GRBs. This leaves dust obscuration of QSOs and 
gravitational lensing of GRBs as the two most plausible 
e xplanations. 

| Wyithe et al.j (2011 hereafter Wll) suggested that, 
being detected in two independent energy bands, GRBs 
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coul d be subjected to a multi-band magnification bias 
(see Wyithe et al. 2003 for details). Consequently, ob- 
served GRBs could be more likely to be strongly grav- 
itationally lensed than QSOs, with a resulting increase 
in absorber number s due to the lensing gala xies. Follow- 
ing this prediction, Rapoport et al.| (|2012) investigated 
the probability distribution ol alignment between GRBs 
with strong Mgll absorption and the galaxy closest to 
their LOS. The most interesting cases to suggest possi- 
ble strong gravitational lensing scenarios were those of 
GRB020405, which showed a second transient ~ 3" from 
the GRB's optical afterglow location, and GRB030429, 
which is aligned only 1.2" away from a massive galaxy 
that is known to have caused strong Mgll absorption. 

As a part of their study of dusty systems along QSO 
LOS JBudzynski fc Hewett]Q2011| hereafter BH11) tested 
the bias in absorber numbers expected due to dust ob- 
scuration by foreground objects towards QSOs and con- 
cluded that while some bias is anticipated, it cannot 
solely account for the discrepancy. From the total sample 
of QSOs they studied, they found that 24%±4% of strong 
Mgll absorbers are not observed because the QSOs are 
obscured and either fall below the typical signal-to-noise 
(S/N) thresholds or are no longer detected. 

In this paper we quantify the significance of the Mgll 
problem when including 3 new LOS, and describe the 
statistical tools we use for this purpose in Section [2] 
In Section [3j we analyze the gravitational lensing bias, 
presenting new observations of the lensing candidate, 
GRB020405; discuss the dust obscuration bias; and cal- 
culate the total bias required to resolve the discrepancy. 
Our conclusions are summarized and discussed in Section 

i 

2. SIGNIFICANCE OF THE MGII PROBLEM 
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To quantify the significance of the Mgll problem we 
require knowledge of the redshift path over which strong 
systems can be observed (also known as path density, 
g(z)) towards the GRBs. The appropriate comparison 
for the number of observed systems towards GRBs is 
the expected number of systems towards QSOs covering 
an equivalent g(z), which is calculated from the QSO 
absorber number density (dn/dz) as 



N, 
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exp,QSO 
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Recently, |Lawther et al.| ( |2Q12| hereafter L12) used 
SDSS DR7 to study QSOs with Mgll absorption lines. 
For each redshift bin they measured how many QSO 
spectra would allow detection of a strong Mgll system 
at that redshift, and the number of absorbers N(z) ob- 
served in that bin. They find that for strong absorbers 
(equivalent width (EW) > lA), the number density per 
co-moving Mpc along a LOS, n(X), can be well described 
as a function of redshift, z, by the following expression: 



cn(X)/H = n exp(z /z), 



(2) 



where no = 0.1 10 ±0.005 absorbers per unit redshift, and 
zo = — 0.11 ± 0.03. The number density, dn/dz is given 
by: 



dn cn(X) 



{l + zf 



(3) 



In 2006, P06 reported that the probability of find- 
ing the number of strong Mgll absorbers observed to- 
wards 14 GRBs in a similar redshift path towards QSOs 
was < 0.1%, using a Monte-Carlo (MC) analysis. In 
2009, V09 analysed all GRBs with high-resolution spec- 
tra taken by UVES on the VLT up to June 2008 (10 ob- 
jects). Their complete sample, which included 26 GRBs 
(16 with available data from the literature, including 
those from the P06 sample), revealed 22 strong Mgll 
systems. They found that for a similar redshift path 
towards QSOs, one would expect 10.41 strong absorb- 
ing systems. However, they incorrectly included Poisson 
statistics to express the error for this number when it 
was calculated from a fit to the absorber number density 
equation (dn/dz) for QSOs. We advocate instead that it 
is best to use a MC approach to reflect the distribution 
of absorbers, such as was undertaken by P06. 

Here, we revise the V09 work, deploying different 
statistical tools, and also include 3 new GRBs with 
high-resolution VLT data (GRB080804, GRB081008 and 
GRB081029). 

2.1. Additional Data 

We reduced and analyzed publicly available data for 
GRB080804, which we acquired from the ESO archive. 
The data were taken using the UVES instrument on the 
VLT, with a total 2.6 hours of integration time. The 
data were acquired with a 1" slit in good seeing (~ Of! 7). 
The observations were reduced using the ESO Gasgano 
pipeline for UVE^j The pipeline performs the standard 

4 http:/ /www. eso.org/projects/dfs/dfs-shared/web/vlt/vlt- 
instrument-pipelines.html 
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Fig. 1. — Redshift path density g(z) for absorbers with EW > 
lA towards the GRBs used in this study. The contribution labelled 
"Vergani 2009" includes only the 10 UVES GRBs, while the values 
from the literature were also included in the full sample in V09. 
The objects we added to the study are shown in yellow. The total 
sample used in the paper is the sum of all the different contributions 
and is outlined with the black curve. 

tasks of bias subtraction, flat-fielding, spectral order lo- 
cation, wavelength calibration and extraction of the spec- 
trum. Candidate Mgll doublets were identified by visual 
inspection of the spectra. The high signal to noise ra- 
tio of the data allowed lines exhibiting EWs as low as 
0.1 A to be detected. All candidate doublets were then 
shifted to the rest-frame to assess if the lines matched 
the Mgll doublet at 2794.4A and 2801. 5A. All equivalent 
widths were measured in IRAF0 using the "splot" task 
(Tody 1993). The total wavelength coverage is 3,600A- 
9,000A with telluric features preventing possible identifi- 
cation of the doublet between 5,598A- 5,670A, 6,850A- 
6,900A, 7, 170 A- 7,350A and 7,520A- 7,670A. The wave- 
length limits were converted into redshift space for inclu- 
sion in the redshift path density. No strong Mgll absorb- 
ing system was identified. 

The redshift path and absorber information for GRBs 
081008 and 081029 were taken fro m the literatur e (ID 'Elia 

respectively). The full 
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dataset is listed in Table [TJ and the corresponding red- 
shift path density is shown in Figure [I] 

2.2. Statistical Analysis 

We calculate the number of strong Mgll absorbers one 
would expect to find towards QSOs for the redshift path 
of our observed GRBs using equation [I] Using the V09 
redshift path density for only the 10 UVES GRBs, which 
included 9 strong absorbing systems, LI 2 predicted that 
4.1^0*7 strong absorption systems would be expected 
from a similar QSO sample. While the relatively small 
errors represent the degree of precision in the mean ex- 
pected value for strong absorbers along QSOs LOS, we 
note that the probability around the mean value is not a 
Gaussian distribution, and therefore that the error can- 
not be used to directly estimate the significance of the 
Mgll problem. 

5 http://iraf.noao.edu 
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Using dn/dz from L12, and the g(z) for our 29 GRBs, 
we find A^| n = 9.6 ± 0.7 for a comparable sample of 
QSO LOS, where the errors are propagated from the un- 
certainties of the fitting parameters for dn/dz (no and 
zq). When using the P06 values for no and zo (values 
found for SDSS DR3), A^q SO = 7.4 ± 0.7 with the 
g(z) provided in V09 only, and 8.3 ± 0.8 for the g(z) in- 
cluded in this work. The difference between the GRB 
and the mean of the QSO populations is reduced when 
using our sample relative to the earlier V09 sample. It 
is further reduced when using the updated L12 fit in- 
stead of the P06 fit to the QSO number density (as the 
expected value for QSO absorbers increases). Thus, as 
more data are available, both for GRBs and QSOs, the 
difference between the expected number of absorbers to- 
wards QSOs and the observed population towards GRBs 
is seen to decline, implying that we are still strongly af- 
fected by small number statistics. 

In order to assess the significance of the current Mgll 
problem we calculate how unlikely it would be to observe 
22 strong Mgll systems along the LOS to 29 GRBs, un- 
der the hypothesis that GRBs and QSOs probe the same 
absorber population . We use the SDSS DR4 catalog of 
|Quider et at] fl2011| ), which includes 44,600 QSOs spec- 
tra and contains ~ 17, 000 measured Mgll doublets. The 
catalog provides information regarding the redshift of the 
QSO and the redshift and EW of the Mgll absorber. We 
conduct a MC simulation in which we randomly select 
QSOs from the sample with similar g(z) to the observed 
GRBs, and count the number of strong absorbers to- 
wards them. We are not constrained to match the precise 
number of source systems, as the absorber rate depends 
only on the redshift path being probed. Therefore, for 
each redshift bin of the g(z) function (see Figure fl| we 
select QSOs at a redshift of a randomly selected (JRB 
(within a 0.1 redshift bin), constraining the GRB sam- 
ple to those at a redshift which is larger than the bin, 
and small enough so the tested redshift bin does not fall 
shortward of the Lyman alpha break in the QSO spec- 
trum. Once an appropriate random QSO is chosen, we 
count the number of absorbers within the redshift bin. 
We continue this procedure until we have covered the 
full redshift path of our g{z) function. Our MC analysis 
includes 10, 000 trials of the full path selection. 

Our analysis reproduces the mean value calculated 
above from the empirical fits of P06 (8.3 absorbers). The 
histogram of the number of absorbers is consistent with a 
Gaussian distribution, for which we measure a standard 
deviation of a = 3.2 absorbers. Our simulation resem- 
bles the results of P06 more than those of L12 as we are 
using the QSO Mgll absorbers catalog from SDSS DR4. 
The SDSS DR4 is the largest publicly available catalog 
which includes all the data required to conduct this anal- 
ysis. P06 found fewer absorbers per unit redshift in the 
DR3 sample than the L12 DR7 study, as seen in Figure 
3 of L12. The observation of 22 absorbers among a pop- 
ulation with a mean of 8.3 and a — 3.2 would be 3.9cr 
from the average value. As expected if the distribution 
is Gaussian, we find no instances of a sample with 22 
or more strong absorbers in our 10,000 MC trails. To 
test the sensitivity of our results to the chosen QSO red- 
shift we repeat the MC analysis and do not constrain the 
QSOs to have similar redshifts as the GRBs. We do not 



find any statistically significant deviations in the results, 
either in the absorber numbers or in the redshift distri- 
bution of the absorbing systems. Moreover, the redshifts 
of the absorbing systems are similar between the QSOs 
and the GRBs. 

3. POTENTIAL BIAS CONTRIBUTIONS 

Having found the difference between the GRB and 
QSO lines of sight to be statistically significant, we now 
explore the two leading explanations for this discrepancy: 
gravitational lensing bias and dust obscuration bias. 

3.1. Gravitational Lensing Bias 

A potential contributor to the Mgll problem is that 
GRBs with observed afterglows are more likely to be 
gravitationally lensed than QSOs (see Wll and refer- 
ences therein). This is attributed to the multi-band mag- 
nification bias, arising from GRBs being detected in two 
bands (7- rays and optical) w ith uncorrelated intrinsic 
fluxes (see iWyithe et al.||2003| for details). If GRBs are 
preferentially lensed compared to QSOs, a higher rate of 
Mgll absorption would be expected towards GRBs from 
the gaseous halos surrounding the lensing galaxies. 

3.1.1. GRB020405 

Following the Wll prediction that, if gravitational 
lensing bias is the explanation for the Mgll problem, then 
10 — 60% of GRBs with strong Mgll absorption should 
have been lensed strongly enough to produce multiple 
images, we studied archival data of the GRBs from the 
V09 sample to look for potential strongl y gravitation- 
ally lensed systems ( Rapoport et al.||2012 ). Two poten- 
tial cases of lensing were identified, one of which was 
GRB020405 (z=0.695), which had a nearby (3") tran- 
sient that could have been a repeating image of the same 
GRB. 

The transient near GRB020405 was first observed by 
Masetti et al. (2003) in an HST image taken ~ 19 days 
after the Swift trigger. They noted that it was not visi- 
ble in F555W but only in F702W and F814W, and had 
completely faded by the last observation in August 2002. 
Figure |2] shows the locations of the GRB, transient and 
host galaxy. Objects 1 and 2 were confirmed, using 
VLT spectroscopy, to be at the s trong absorbe r redshift 
(z = 0.472). Using GRAVLENS ( |Keeton|20Ql] ) to model 
the system, we found that if objects 3-6 (or even only 4-6) 
are part of a group at the absorbers' redshift, the obser- 
vati ons could be explaine d as a multiply imaged GRB 



(see Rapoport et al. 



2012 



for details and Figure [2] right 
panel tor a potential model). In order to test this sce- 
nario we initiated a detailed analysis of the surrounding 
galaxies. 

We obtained spectra for several of the gal axies around 
GRB 020405 using the GMOS instrument flHook et al~ 
2004) on Gemini-South (PI Rapoport, Program ID GS- 
^UT2A-Q-9). The multi-slit data were taken on UT 03 
April 2012 and consisted of eight 120s exposures with 
the R150 grating. The 1" slit widths provided a spectral 
resolution of R ~ 315, and the data were binned 4x4 to 
yield a spatial scale of ~ 0.3" pixel -1 and a dispersion 
of 13.6A pixel- 1 . To improve the removal of sky and 
fringing features, the mask was cut with two sets of slits, 
offset by 85", and the field was dithered between the two 
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positions. The data were reduced using the standard 
Gemini/IRAF packageer] for flat fielding and wavelength 
solution. No standard flux correction was applied as the 
only information needed from the observations were red- 
shifts. 
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Fig. 2. Left: HST WFPC2/F702W field of GRB020405. The 
GRB is clearly visible and the complex host and second transient 
are indicated. Objects 1 and 2 were found to be at the strong ab- 
sorber's redshift of z=0.472. Our observations indicate that object 

3 is a field star, galaxy 4 is a galaxy at z = 0.485 ± 0.002, object 5 
has unknown redshift and galaxy 6 is at z = 0.484 ±0.002. The slits 
used in the Gemini/GMOS observations are superimposed on the 
image. Right: LENSMODEL solution suggested in Rapoport et 
al. (2012), which assumed objects 1-6 were at the strong absorber 
redshift of z=0.472. With the newly measured redshifts of objects 

4 and 6, we are not able to identify any gravitational lensing model 
to explain the transient as another image of the GRB. 



Comparing our reduced spectra to a library of tem- 
plates using the RUNZ 2dFGRS redshift code (R. Sharp 
2012, private communication, see Figure [3|, we were 
able to determine that objects 4 and 6 are galaxies at 
z = 0.485 ± 0.002 and z = 0.484 ± 0.002 respectively, 
where the errors on the redshifts are based on a centroid 
fit to the strongest line. Object 3 is consistent with being 
a foreground star. The spectra for object 5 were not of 
sufficient signal-to-noise to determine the type or redshift 
of the object. 

As the galaxies at z = 0.472 and z = 0.484 are too 
distant to be from the same halo (At; = 4800 km/s), but 
objects 4 and 6 are at a similar redshift, we modelled 
the field using two different groups (objects 1 and 2, and 
objects 4 and 6). In order to simplify the mathemat- 
ics we approximated the two groups (each as a singular 
isothermal ellipse with shear) to reside in one plane at 
the strong absorber redshift (z = 0.472). No model was 
found which could explain the observables using these 
galaxy redshifts and allowing plausible masses for the 
galaxies. 

Therefore, we conclude that the transient near 
GRB020405 was not another image of that GRB, leav- 
ing one potential candidate as a strongly gravitationally 
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Fig. 3. — GMOS spectra of galaxies 4 (lower panel) and 6 (upper 
panel), nearby galaxies to GRB020405 (see Figureplfor alignment). 
The strongest emission lines are identified and firm galaxy 4 to be 
at z = 0.485 ± 0.002 and galaxy 6 at z = 0.484 ± 0.002. 



lensed GRB from the sample studied in our previous 
work, GRB030429. This GRB was detected 1.2" away 
from a strong Mgll absorbing spiral galaxy. Using pho- 
tometric observations, a spectral energy distribution fit- 
ting technique implied Mb = —21.1 ± 0.1 for the ab- 
sorbing galaxy. The Tully-Fisher relation for such an 
intrinsically bright galaxy suggests a velocity dispersion 
of 160 ± 65 km/s, where 200 km/s is r equired for pro- 
ducing a second image for this GRB (see |Rapoport et aL] 
[2012] for further details). Thus, GRB030429 is likely to 
be magnified at some level, possibly strongly. Verifying 
this scenario would require measuring the velocity dis- 
persion of the galaxy, w hich with R = 22.70 ± 0.12 mag 
( Jakobsson et al.|[2004 ) and at a redshift of 0.841 would 
be observationally very expensive. 

In our previous analysis of the existing HST imaging 
data for 11 GRBs with Mgll absorption from the V09 
sample, we ruled out strong gravitational lensing in 6 
cases. Now having determined that GRB020405 was not 
strongly lensed, we can place an upper (2a) limit on the 
fraction of strong absorbers which lead to multiply im- 
aged GRB of < 35%, though this is a very conservative 
limit as the remaining systems show no evidence of lens- 
ing. Assuming only the most likely case of GRB030429 
is a viable strongly gravitationally lensed candidate, the 
lensing fraction is reduced to < 13% (2a). 
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3.1.2. Statistical Analysis 

To test the effect of strong gravitational lensing on 
the Mgll absorber statistics, we account for the frac- 
tion of strongly gravitationally lensed GRBs (F\ ens ) by 
multiplying the number of absorbers found in the MC 
trials by 1/(1-Fi ens ). This approach assumes that all 
strongly gravitationally lensed objects would show strong 
Mgll absorption. We find the probability of finding 22 
absorbers for a lensing fraction of 0.05 (assuming the 
Mgll absorber towards GRB030429 was the only one 
out of the 22 systems that caused strong lensing) to 
be 0.02 %. Moreover , even if all 4 remaining candidates 
from Rapoport et al. (2012) were strongly gravitationally 
lensed, the probability ol observing 22 absorbers would 
only increase to 0.7%. This indicates lensing alone can- 
not account for the Mgll problem. 

3.2. Dust Obscuration Bias 

The Mgll problem could also arise by flux-limited QSO 
samples selecting against absorbed objects with high 
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Mgll EW. Among a QSO sample having the same red- 
shift path density as GRBs in the V09 sample, BH11 
showed that ~ 20% (see Figure 19 in BH11) of sources 
with strong absorbing systems would have been obscured 
by dust and so either not detected or below the S/N 
threshold (for EW > lA^j Thus, the QSO sample is 
likely missing 20% of the strong absorbing systems. To 
correct for this dust obscuration, we multiply the pre- 
dicted absorber rate from our MC analysis by a factor 
of 1/(1 — 0.2). This correction results in a mean of 10.8 
absorbers, with a probability of finding 22 absorbers of 
0.4%. 
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Fig. 4. — Redshift s of the QSO s , GRB s and their strong absorb- 
ing systems. The |Quider et al.| ( |2011| ) QSO sample used in the 
MC analysis is represented by grey points and their corresponding 
normalised histograms. The black triangles and the black stars 
mark the GRBs with absorbers from the full sample used in this 
paper, where the stars differentiate those GRBs which would have 
been excluded in the dust obscuration analysis of BH11. The lower 
black line represents the ZMgll — zqso boundary and the upper 
black line is the limit due to Lyman alpha absorption. 



3.3. Total Bias 

As neither dust obscuration nor gravitational lensing 
can by themselves explain the Mgll problem, we repeat 
the MC simulation and investigate the total bias required 
to resolve the difference. We do this by imposing a bias 
on the QSOs by multiplying the number of absorbers 
by 1/(1 — FBias)- We then increment the bias fraction 
(^Bias) until it can account for the observed excess to- 
wards GRBs. We find that for a total bias fraction of 

7 At higher redshift one would expect the QSOs to appear fainter 
and be more susceptible to dust obscuration. As the added g(z) in 
this work, relative to that in V09, is flat (see Figure!!]), the redshift 
path could be thought as being more concentratecfat higher red- 
shifts than it is in the V09 sample. Therefore, by using the value 
found by BH11 for the percentage of obscured systems, we are more 
likely to be underestimating the number of obscured QSOs than 
vice- versa. 



~ 30% (45%) the significance of the problem is reduced 
to 2cr (la) (see Figure [5]). This formalism assumes that 
the majority of the bias arises from the QSOs, which is 
consistent with the ~ 20% bias estimated by BH11. 

BH11 noted that because their sample is flux limited 
they are likely to be underestimating the fraction of ab- 
sorbers missing from the SDSS DR7, with systems hav- 
ing EW > 5. OA being completely removed. Moreover, in 
their analysis they only use QSOs up to a redshift of 3.5 
and limit the number of absorbers per LOS to two. Fig- 
ure [4] shows the QSO, GRB and absorber redshifts used 
in our MC analysis. To demonstrate the effect BHll's 
QSOs selection would have on the GRB sample, we flag 
the GRBs that would have been omitted with the star- 
like markers. The large fraction (60%) of absorbers along 
the LOS to GRBs which would have statistically gone 
through more dusty systems than the QSOs that were 
examined by BH11 illustrates that the dust obscuration 
bias could be underestimated. Therefore, while strong 
gravitational lensing might be playing a minor role, the 
likelihood that the dust obscuration bias has been under- 
estimated implies that this effect can more easily resolve 
the Mgll problem. 
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Fig. 5. — Probability of finding 22 strong absorbing systems along 
the LOS to QSOs having a similar redshift path density to that of 
the GRB sample for different fractions of total bias. The number 
of absorbers found along the QSOs LOS in the MC is multiplied 
by l/(l-i ? Bias) to account for the missing fraction before being 
compared to the number of absorbers found towards GRBs. 



4. DISCUSSION 

We have presented an updated analysis of the GRB 
strong Mgll absorber problem, including additional GRB 
lines of sight. We outline a statistical method with which 
the Mgll problem can be addressed, using a MC tech- 
nique to estimate the number of absorbers expected to- 
wards QSOs for an observed redshift path towards GRBs. 
We find that the observed discrepancy between the num- 
ber of strong Mgll absorbers in GRBs and QSOs has a 
level of 4<r significance, with < 0.01% chance of finding 
the observed number of absorbers along GRB lines of 
sight in the MC results for the QSOs. 

Gravitational lensing bias for GRBs can contribute to 
the discrepancy, and we examined the best lensing can- 
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TABLE 1 
GRB SAMPLE 



GRB 


^GRB 1 


Az 2 




Reference 4 


yyizio 


1.022 


0.636 


0.770 
0.803 


1 


UUU9ZD 


2.038 


1.392 




2 


010222 


1.477 


1.022 


0.927 
1.156 


3 


011211 


2.142 


1.566 




1 


020405 


0.695 


0.312 


0.472 


4 


020813 


1.255 


0.866 


1.224 


5 


U21UU4 


2.3295 


1.756 


1.3800 
1.6026 


6 


nonooc 

UoUzzo 


1.986 


1.590 




7 


U3U323 


3.372 


0.822 




8 
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2 Redshift path length for strong Mgii absorbers 
^ Strong absorbe r redshift 

4 References 1; |Vreeswijk et al.| (|2006[ 2: |Castro et"aLl ( |2003[ 3: 
Mirabal et al.| ( |2Q02| ) 4: |Masetti et al.|fl2Q03| ) 5; |Barth et al.|(|2QQ3|) 6 : 
Vergani etTTH2009| 7: |Klose et al.| ^ 2 004^8: [yTeeswijk et al.|fl2004|) 
|Maiorano et a l. (2006) 10: Jakobsson et al. (2 004) 11 : |Berger et"af] 
2006 ) 127 ]Prochter et al. ( 2006b )" l3; |Chen et al.|(|2009| ) 1 4: |Hao et all 
2007| 15: Thone et al.|fl2008|) 1 6: |Perley et al.| ( |2008| ) 17: ID'Elia et al"] 
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didate of [Rapoport et al.| ( |2Q12] ) ? GRB020405, which has 
a known foreground absorber and a second nearby tran- 
sient. We find that the nearby galaxies, which would 
need to be at the absorber's redshift to allow possible 
strong gravitational lensing by a galaxy group, are not at 
the appropriate redshift. With no feasible lensing model 
found for the GRB and its nearby transient, we con- 
clude that GRB020405 was not affected by strong lens- 
ing. However, there remains the possibility that gravita- 
tional lensing is playing a role, as evidenced by the very 
close alignment of GRB030429 to a massive foreground 



galaxy. If this GRB is significantly lensed, then the Mgii 
discrepancy is slightly reduced, with the probability of 
finding the number of strong absorbers along the GRBs 
inc reasing to 0.02% . If 4 G RBs from the sample studied 
by | Rapoport et al. ( |2012| are strongly lensed, this im- 
plies a lensing traction ol 20%, which would increase the 
probability of finding the observed number of absorbers 
in a comparable QSO sample to only 0.7%. Thus we 
conclude that lensing is very unlikely to solve the Mgii 
problem on its own. 

The significance of the problem is also reduced after 
accounting for dust obscuration bias, with the probabil- 
ity of finding the observed number of absorbers along the 
LOS to the GRBs increasing to 0.4%. When including 
both dust and gravitational lensing biases, we calculate 
a probability of 0.9%, which remains statistically signifi- 
cant. In order to reduce the problem to below a 2<j level, 
a total bias of ~ 30% is required. We suggest this could 
be largely satisfied by the dust obscuration bias towards 
QSOs currently being underestimated. 

Following submission of our paper, Cucchi ara et al.| 
fl2012| ) increased the sample of GRBs to 118 LOS by in- 
cluding GRBs with low-resolution spectra. This study 
has found the rate of Mgii absorbers to GRBs to be lower 
than previous studies, reducing the overall discrepancy 
between the absorber rate to GRBs and QSOs to < 90% 
confidence level. If additionally, one were to include the 
dust obscuration bias described above, we expect the rate 
of absorbers between GRBs and QSOs to be completely 
consistent. 
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